Identifying Patterns on Cosmic Ray Maps with Wavelets on the Sphere 
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The deflection of ultra-high energy cosmic rays depends on the shape of the injection spectrum of the source 
and the pervasive cosmic magnetic fields. In this work it is applied the wavelet transform on the sphere to 
search for energy ordered filamentary structures arisen from magnetic bending. These structures, the so-called 
multiplets, can bring relevant information concerning the intervening magnetic fields. 
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I. INTRODUCTION 

Ultra-High Energy Cosmic Rays (UHECRs) are particles with 
energy above 10 18 eV reaching the Earth. Almost a century 
after the discovery of cosmic radiation, the origin, chemical 
composition and mechanisms of acceleration and propagation 
of the UHECRs remain a mystery (TJ. 

The way their arrival directions are represented is by associ- 
ating each event with its corresponding latitude and longitude 
in the celestial sphere. For charged particles, the pervasive 
cosmic magnetic fields play an important role, affecting the 
distribution of arrival directions of UHECRs. 

Several mechanisms of acceleration of UHECR in different 
astrophysical sites result on a power law differential energy 
spectrum[3|, i. e., 
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where a is the spectral index. For a single source emitting 
UHECRs of energy E according to this spectrum, the expected 
angular deflection is inversely proportional to the energy. The- 
refore, particles with different energies will deflect differently 
and the greater the energy of the particle the smaller the an- 
gular deflection of the particle (and vice-versa). So, if locally, 
the typical scale of the intervening magnetic fields (coherence 
length) is smaller or of the same order of magnitude of the 
Larmor radius of the cosmic ray, filamentary structures in sky 
maps are expected to be formed. The cosmic ray events that 
compose this structure are ordered by energy, creating the so- 
called multiplets of UHECRs. For the galactic magnetic field 
the expected angular deflection (8) is given by[2] 



5 = 8.1°40- 



B 



dr 



3 kpc 2 kpc 



(2) 



where L is the distance to the source, B is the magnetic field, 
Z is the atomic number of the nuclei and E its energy in EeV. 

II. THE METHOD 

In this work it is proposed a new method (for another method 
refer to (4]) to identify multiplets in maps containing arrival 
directions of UHECRs. This method is based on the spherical 



wavelet transform! 5 1 an d consists on the correlation of a sig- 
nal with rotated versions of a given pattern, both defined on a 
spherical manifold. Let a\ m be the coefficients of the spherical 
harmonic expansion of the signal, and b[ n be the same expan- 
sion for the sought rotated pattern (in this case, a wavelet). 
The correlation coefficient is 
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where the overline indicates the complex conjugate, D l , de- 
signates Wigner-D functions, B is the band limit of the imple- 
mented algorithm, and (a,/3, y) are the Euler angles associa- 
ted to the rotation of the wavelet. 

Applying the spherical wavelet transform using a code cal- 
led SWAT (Spherical Wavelet Analysis Tool)||6l to the dis- 
tribution of arrival directions of UHECRs, if the correlation 
coefficient is high, this indicates a similarity between the pat- 
tern (wavelet) used and the signal. The Euler angles provide 
information regarding the coordinates (latitude and longitude) 
where there is a maximal correlation coefficient, and the ori- 
entation of the pattern. 

Even though the method is able to identify filamentary 
structures, it does not take into account the energy ordering of 
the events for the analysis. This is done by choosing a rectan- 
gle of dimensions ~ 12° x 1.5° (typical dimensions of a mul- 
tiplet) with the obtained orientation, and calculating the Pear- 
son coefficient for 8 x E^ 1 of the events within this rectangle. 
If this coefficient is high, it indicates an energy ordering and 
the filamentary structure is a possible multiplet candidate. 

III. APPLICATION TO SIMULATIONS 

Events of UHECRs were simulated using the code CRT|[7), 
which numerically solves Lorentz force equations of motion 
for charged particles in a magnetic field using fifth order adap- 
tative Runge-Kutta routines. Fifty proton events were simula- 
ted from a hypothetical source lying at a distance of 4 kpc 
from Earth, with galactic coordinates (l,b) = (90°,— 45°). 
The source was assumed to have an injection spectrum with 
a =2.7 and to emit particles within a perfectly collimated jet. 

The magnetic field models used in the simulation are the 
ones proposed by Harari, Mollerach and Roulet[8|, namely 
the ASS-S, ASS-A, BSS-S and BSS-A models, were ASS 
stands for AxiSymmetric Spiral (even under the transforma- 
tion 9 —> 9 + K, where 9 is the azimuth angle) and BSS for 
BiSymmetric Spiral (odd with respect to the mentioned trans- 
formation). The sufixes -A and -S indicate the parity of the 
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transformation above and below the galactic disk, where -S is 
symmetric (even) and -A is antisymmetric (odd) with respect 
to the z — > z transformation. In this class of models the mag- 
netic field in the disk of the galaxy can be decomposed in a 
radial and an azimuthal part, depending on the pitch angle (p) 
of the logarithmic spiral that models the magnetic field of the 
Milky Way. The intensity of the magnetic field in the disk is 
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where ro = 10.55 kpc, r\ = 2 kpc, p = —10° and m is 1 for 
BSS models and 2 for ASS models. The z component of the 
magnetic field is 
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Figura 1: Skymap containing the arrival direction of the 
simulated events. The size of the circles are proportional to 
the energy of the event. 



where £(z) = tanh(z/z2) for -A models and 1 for -S models, 
zo = 4 kpc, z\ — 0.3 kpc and zi — 20 pc. 

By applying the method to these specific simulations, using 
the wavelet transform with suitable parameters (compatible 
with the expected dimensions of a multiplet), the multiplet 
was correctly identified and the orientation for different mo- 
dels was calculated and is shown in table|I] Figure[TJshows the 
simulated events for the simulations. 

Tabela I: Orientation (%) for the simulated source in the class 
of models of Harari, Mollerach and Roulet, for the different 
symmetries. 

Model x(°) 
ASS-S 100.3 
ASS-A 108.4 
BSS-S 65.7 
BSS-A 68.1 



IV. CONCLUSIONS AND PERSPECTIVES 

The orientation of the multiplet, i.e., the angle it forms with 
the galactic plane, depends on the coordinates of the source 
and the magnetic field model adopted. Therefore, it is possi- 
ble to set limits on models of the galactic magnetic field by 
analysing the orientation of multiplets. When the method was 
applied to simulated data, the multiplets were correctly iden- 
tified and it was noticed a relation between the orientation of 
the multiplet and the adopted model. If this method is ap- 
plied to data collected by the Pierre Auger Observatory or any 
other cosmic ray experiment and is able to successfully iden- 
tify multiplets, a new tool to probe the galactic magnetic field 
will be available. 
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